The effectiveness of HYDRAULIC FRACTURING fluid injection is influenced by numerous factors, including pre-existing discontinuities such as discrete fracture networks (DFNs). Among the geometric characteristics of DFNs, fracture density is a critical factor. In deep reservoirs, which often consist of hot dry rock (HDR), thermal conduction through the rock and fluid, as well as advection and convective heat transfer within the fluid, can significantly impact fluid–rock interactions. This study examines the influence of DFN density on HYDRAULIC fracture (HF) propagation in HDR, with a particular focus on the thermo-hydro-mechanical (THM) behavior of HDR using the combined finite-discrete element method (FDEM). Key controlling factors, such as flow rate, fluid kinematic viscosity, in-situ stress magnitude, pre-existing fracture aperture, and working fluid temperature, are analyzed. The findings highlight the significant role of DFN density in determining the pattern and extent of HF propagation under varying conditions. Additionally, the interaction between the working fluid and DFNs is shown to vary considerably with changes in these controlling factors. However, the study reveals that variations in DFN density or the values of the controlling factors have minimal impact on the temperature field. This is attributed to the rapid heat exchange between the cold fluid and the HDR, which quickly raises the fluid temperature, resulting in negligible temperature variations.